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1 Introduction 



o 

Whenever bending of light by gravitational fields is detected in astronomical 
. observations it is generally very weak. Even in the spectacular formation of 

giant arcs, typical of the so-called strong lensing, photons are deflected by no 
more than a few arcseconds. This leads to the conclusion that the first order 
Einstein formula for the deflection angle is sufficient to explain all observed 
phenomenology. On the other hand, the existence of very compact objects, 
such as neutron stars or black holes, is now well-established through many 
independent astrophysical observations. In the neighbourhood of such objects, 
^jTj' electromagnetic radiation is generated and travels through very strong gravita- 

tional fields. In such extreme cases, the calculation of the deflection of photons 
needs to be pushed beyond the first order approximation on which the Ein- 
stein formula relies. It must be said that no direct observation of gravitational 
lensing by black holes or other compact objects has ever been performed up 
to date. Yet, some peculiarities of the emission spectra of black holes, such as 
the shape of the K a line of the iron, have been interpreted as the final effect 
of the strong deflection of photons emitted by the accretion disk. 

The aim of this article is to review the theoretical aspects of gravitational 
lensing by black holes, and discuss the perspectives for realistic observations. 
We will first treat lensing by spherically symmetric black holes, in which the 
formation of infinite sequences of higher order images emerges in the clearest 
way. We will then consider the effects of the spin of the black hole, with 
the formation of giant higher order caustics and multiple images. Finally, we 



will consider the perspectives for observations of black hole lensing, from the 
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detection of secondary images of stellar sources to the interpretation of iron 
K-a lines and direct imaging of the shadow of the black hole. 



2 Spherically Symmetric Black Holes 

The first presentation of the exact calculation of light deflection by a com- 
pact body using the Schwarzschild metric dates back to the work by Charles 
Darwiifl in 1959 pQ. 

A generic static spherically symmetric metric can be put in the form 

ds 2 = A{r)dt 2 - B(r)dr 2 - C{r){dd 2 + sin 2 M<j) 2 ). (1) 

In particular, the Schwarzschild metric has A(r) = l — 2M/r, B(r) = [j4(r)] _1 , 
C(r) — r 2 . Suppose that a photon comes from infinite distance, grazes the 
black hole at a minimum distance r m and goes away to infinity. The deflection 
angle, defined as the angle between the asymptotic incoming and outgoing 
trajectories, can be easily derived by the analysis of the geodesies equations 

fYl I ' 

SW *, (2) 




C(r) [C(r)/A(r)-u 



where u is the impact parameter, defined as the distance between the black 
hole and each of the asymptotic photon trajectories, related to the closest 
approach distance r m by 

(3) 



2 _ C(r m ) 



A(r m ) ' 

The integral © was solved by Einstein assuming that r > r m 3> 2M. In 
this limit it reduces to the standard Einstein result for the Weak Deflection 
Limit (WDL) 

awDL =4M/u, (4) 

with u ~ r m . Darwin, instead, made the integration exactly in the case of a 
Schwarzschild metric. His deflection angle is 



a = -n + 4\/r m /sF(ip,m), (5) 
where F{<j>, m) is the elliptic integral of the first kinc0, and 



s = y/{r m -2M)(r m + 6M) (6) 
m = (s-r m + 6M)/2s (7) 



ip = arcsin ^2s/{3r m - 6M + s). (8) 



1 Charles Darwin was the grandson of the famous evolutionist with the same name, author 
of "The Origin of the Species" . 

2 In our notations F(tp,m) = f . d ® 

q y 1-m sin 2 9 
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Fig. 1 Comparison between the exact deflection angle for photons in a Schwarzschild 
metric and the two limits discussed in the text. The impact parameter is in units of the 
Schwarzschild radius. 

Eq. ([5]) gives the deflection angle as a function of the unobservable mini- 
mum distance r m . We can better express it in terms of the impact parameter 
u (which is coordinate- independent), using Eq. ([3]). The exact deflection angle 
is plotted in Figure [TJ along with its approximations in the WDL ((3]), and the 
Strong Deflection Limit (SDL) to be introduced below ([§]). 

Darwin's deflection angle and the Einstein approximation coincide at very 
large impact parameters. At small impact parameters, however, the salient 
feature of the exact deflection angle is the divergence at a finite value u = u = 
3-\/3M, corresponding to a closest approach r rn = f = 3M via Eq. ([3]). Light 
rays with impact parameters closer and closer to u experience larger and larger 
deflections even exceeding 2n. This simply means that a light ray deflected 
with a in the interval [2m:, 2(n + l)ir] performs n complete loops around the 
black hole before leaving it. A light ray with impact parameter exactly equal 
to u would perform an infinite number of loops approaching r = f indefinitely. 
Nonetheless, circular orbits of light rays with radius r — f are admitted by 
geodesies equations, but are unstable against small perturbations, which would 
finally drive photons into the black hole or toward spatial infinity, f is also 
called the radius of the photon sphere (for a general definition see Ref. [3]). 
Finally, photons with impact parameters smaller than u are just captured by 
the black hole and fall into the horizon. 

A good approximation for the deflection angle near the divergence is ob- 
tained by expanding the elliptic integral for r m ~ f. We get 



a S DL =-\og(u/u- 1) + log 216(7-4^/3) 



(9) 



which is the Strong Deflection Limit approximation plotted in Fig. [T] 
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Darwin himself intuited the possibility that a source lensed by a compact 
"star" with radius smaller than r would be replicated in two infinite sequences 
of images on each side of the lent[f|. In fact, thanks to the divergence of a, we 
can always find light rays with an impact parameter u close enough to u so as 
to perform n loops around the black hole (in either way) and finally reach the 
observer. These rt-loop images or higher order images (called ghosts by Darwin 
and sometimes mirages or relativistic images) are missed in the WDL but can 
be well-described using the SDL approximation for the deflection angle. 

In the following subsections we will discuss the black hole lens equation 
and derive the general features of the images. Finally we will review the main 
approximation schemes (WDL, SDL and non-perturbative). 



2.1 The Black Hole Lens equation 

A lens equation is a relation between the geometry (positions of source, ob- 
server and lens in a given coordinate system) and the position of the images 
in the observer's sky. Solving this equation, we can get the positions of the 
images as functions of the geometric parameters of the system. 

The position of the observer and the final direction of the photon hitting 
the observer can be considered as initial conditions for the geodesies equations. 
Tracing them backward in the affine parameter, we can write down a formally 
exact lens equation [S]. This approach has been applied to the Schwarzschild 
lens [9] and then generalized to all spherically symmetric spacetimes [10] (see 
also [IT]). 

Let us consider a source at coordinates^ (rs, (f>s) and an observer at (ro, <fro), 
with ro > rs- In typical lensing situations, the photon is emitted from the 
source, approaches the black hole to a minimum distance r m and then escapes 
toward the observer. Therefore, we can divide the trajectory into an approach 
phase (from rs to r m withal f < 0) and an escape phase (from r rn to ro with 
f > 0). By simple integration of the geodesies equations, we can calculate the 
total azimuthal shift experienced by the photon as 



(u,ro,rs) 



BW f< fr. (10) 



C(r) [C(r)/A(r) - u 2 } 



The two integrals separately cover the approach and escape phases of the 
trajectory, u is the asymptotic impact parameter, related to r m by Eq. ([3]). 

However, we must also consider the possibility that the photon directly 
goes from the source to the observer always moving away from the black hole 



3 The existence of higher order images and /or the SDL expansion of Eq. J2J have been in- 
dependently re-discovered several times 4 5 6 7 , demonstrating how easy is to lose memory 
of old papers. 

4 The spherical symmetry ensures that the motion takes place on a fixed plane. Choosing 
our coordinate system so that "So = = "7 2 we get 1? = n/2 during the whole trajectory. 

5 The dot denotes derivative with respect to the affine parameter. 
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Fig. 2 Generic lensing geometry by a static spherically symmetric body L, illustrating the 
quantities defined in the text. The source is in S and the observer in O. 

with r > (e.g. when the source is located between the lens and the observer). 
In this case, there is only an escape phase and the azimuthal shift is 



Note that the two expressions coincide when rs — ► r m , i.e. when the photon 
is emitted in a tangential direction r = 0. 

The lens equation can be written by simply imposing that the azimuthal 
shift equals the difference in azimuth between source and observer modulo 2tt. 



where <P = <P- if Acf> > <P±(u,ro,r m ) and <P = <P+ otherwise. The integer 
number n takes into account the fact that is a periodic coordinate. 
The angle 9 at which the observer detects the photon is [T2] 



Once we fix the radial coordinates of source rs and observer ro , and their 
relative azimuthal position A<p, by inverting the lens equation (|12p we can get 
u and thus the positions of the images 9 through Eq. (|13l) . 

An interesting limit is obtained when source and observer are very far from 
the lens with respect to all other distance scales (it and M). Taking only the 
smallest powers in u/ro and u/rs, we get the Ohanian lens equation [5] 




(11) 



r ,r s ) = A(j) = <j) Q ~ c^s + 2mr, 



(12) 




(13) 



a- 6 -6 S 



7 = Alp — 7T, 



(14) 
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Fig. 3 Position of the images as a function of the source position A<f> for rg = 2000M (a) 
and rg = 2.6M (b). In both cases tq = oo. e = (0/8 — 1) is the relative distance of the 
image from the shadow border. 



with sin# ~ u/ro, sin 83 = u/rs and a being the deflection angle between 
the asymptotic incoming and outgoing trajectories (a — ^(u, 00,00) — it, see 
Fig. The Ohanian lens equation is an appealing simplification of Eq. (|12p . 
in that it gets rid of the two parameters ro and r$ ■ While the limit ro — 00 
is justified from the fact that the Earth is very far from any known candidate 
black hole, the limit rj » h leaves out the very interesting case in which 
the source is part of the environment of the black hole. A discussion on the 
errors committed when using the Ohanian or other approximate lens equations 
appeared in the literature is contained in Ref. [T5] . 



2.2 Position of the images 

Armed with our exact lens equation ((T2|) . we can now find the images for a 
generic source position varying from A(f> = (source in front of the black 
hole) to Atf> = 7r (source behind the black hole). Only images on the same 
side of the source (direct images) can be obtained by Eq. ([121 , since we have 
assumed that <f> is always positive. However, thanks to the spherical symmetry, 
the position of the images on the opposite side of the source (indirect images) 
can be obtained by solving Eq. (|12l) with A<f> — > 2tt — A<\>. Thanks to the 
2mr on the right hand side of Eq. (|T2"j) . we have one image in each interval 
(m — < <P < mir for any integer m. We will refer to this integer as the 
order of the image. The primary image will thus have < <P < it, the secondary 
image it < <P < 2ir and so on for higher order images. Odd order images are 
direct and even order images are indirect. 

In Fig. [21 we plot the position of the images as a function of Acf). The 
left panel, represents the case of a distant source {r$ — 2000M), while the 
right panel considers a source very close to the black hole between the horizon 
r = 2M and the photon sphere r = f — 3M. The position of the images 
is expressed in terms of e = (9/8 — 1), where 8 = u/ro is called the shadow 
border for reasons that will be clear in a while. From top to bottom, the curves 
represent the primary image (n = 0), the secondary image (A<f> — > 2tt — Acf>, 
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n = 0), the third order image (n = 1), the fourth order image (A<f> —> 2tt — Acf), 
n — 1), the fifth and sixth order images (n = 2). 

The familiar gravitational lensing geometry with the source behind the 
black hole is recovered at A<j> = tt. In case of perfect alignment, we expect the 
primary and secondary images to merge in the classical Einstein ring. Indeed, 
in Fig. la we see that they tend to the same value of e, that is to the same 
angle 9 (but on opposite sides). As soon as the source moves out of perfect 
alignment (Acf) < tt), the primary image moves farther from the black hole 
while the secondary image gets closer. For a source in quadrature (A<f> ~ 7r/2) 
the light bending is practically negligible for the direct image, so 9 reaches the 
maximum value rgfro and then decreases to zero when A<j) = (source in 
front of the black hole). On the other hand, the secondary image gets closer 
and closer to the black hole but always stays outside the circle of radius 9. 

Coming to higher order images, we see that they also form Einstein rings 
when A<j> = tt, but with radius just very slightly larger than 9. As A<p decreases, 
the direct images move farther and the indirect get closer to the shadow border. 
Without the logarithmic scale, all images would appear packed together at 
the shadow border 9. Finally, when the source is in front of the black hole 
{A(j> = 0), every indirect image merges with the direct image of the next 
order to form an Einstein ring. Therefore, black holes surprisingly behave in 
a similar way with sources in front of the black hole and behind the black 
hole, since we can have an infinite sequence of concentric Einstein rings in 
both cases. The geometry with the source in front of the black hole is also 
known as retro- lensing 14,15]. It is important to stress that in no case images 
of sources outside the photon sphere can appear inside the circle of radius 9, 
which justifies its name as shadow of the black hole. 

Now let's move to Fig. which represents the same diagram for a source 
inside the photon sphere but outside the horizon. In this case, e — 9/9 — lis 
always negative, which means that the images form inside the shadow. The 
distance order of the images from the black hole is now reversed, with lower 
order images being closer to the black hole and higher order images appearing 
closer and closer to the border of the shadow. Apart from this major difference, 
the diagram of Fig. is quite similar to with the formation of Einstein 
rings at A<j) = 0, tt. Just note that at A<f> = the direct image is again in 
9 = 0, corresponding to e = — 1. 

The dependence of the images on the source distance is illustrated in detail 
in Fig. |U This figure shows how the radius of the Einstein rings (in standard 
lensing configuration A(f> = tt) depends on the source position. In Fig. [4^ 
we have the primary Einstein ring and in Fig. |4Jd we have the second order 
Einstein ring. In both cases the rings enter the shadow when rs < 3M. At 
very large rs, however, the primary Einstein ring grows according to the WDL 
formula 9e — ^/4:Mrs/ro(ro + rs), whereas the higher order Einstein rings 
stay close to the shadow border. 
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Fig. 4 Position of the Einstein ring vs source distance, (a) Primary ring, (b) Second order 
ring. The horizontal line represents the shadow border u = 3v3M. 



2.3 Magnification 

The magnification of the images only makes sense for sources in the asymp- 
totically flat region r$ 3> M . In fact, in order to compare angular elements of 
image and source, we need to project the source element to a Minkowski space. 
This process is coordinate-dependent unless the source is in the asymptotic 
region. Furthermore, the intensity of the image elements is also red-shifted if 
the source is deep in the gravitational field of the lens as a consequence of the 
conservation of the number of photons Io/vq — Is/ V %- 

Let us consider a source far from the black hole. An angular area element of 
the image has size du/ro in the radial direction and udrj/ro in the tangential 
direction, where 77 is an angular coordinate around the axis lens-observer. The 
corresponding tangential displacement in the source surface is r$ sin A<p drj. 
The tangential magnification is therefore 

u l r ° 

= . A ,,, — , (15) 

r s sin A(p/dos 

where 



dos = y r o + r s ~ 2r ors cos A</> (16) 

is the (Euclidean) distance between source and observe:^. 

For the radial displacement the situation is more subtle. A shift in the 
impact parameter du selects a geodesic close to the original one. By differen- 
tiation of the lens equation (fl"2"|) . we get 

—du + g^dr s = -d(f> s - (17) 

This equation approximates the perturbed geodesic in a neighbourhood of 
the source (rs,<t>s)- The point in which this geodesic will intercept a section 
orthogonal to the emission direction x l e is (rs, </>g) + dy l with 



dy i = (dr s ,dcj)s) = (u/r s ,T^ - u 2 AIAs^) dR, 



(18) 



6 The Euclidean distance can be safely defined for sources and observers in the asymptotic 
region r$,ro S> M, otherwise it loses meaning. 
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Fig. 5 Magnification of the images as a function of the azimuthal shift Atj> for the 
Schwarzschild lens. 



which satisfied r\ijx\dyi = and r\ijdy l dyi = dR 2 . The distance dR spans a 
source element. We can thus define the angular source element in absence of 
gravitational lensing as dR/dos- The double sign in d(j)s in Eq. (fl%|) follows 
the sign of <P± that applies to the image under examination. 
Finally, the radial magnification is 

dujro 

/ir = dRjdoS' (19) 

where du and dR are related by Eqs. (jT7)) and (fT5|) . The total magnification 

is fl — fl r flf 

When ii < rj and the photon is emitted with f < 0, our definition of 
magnification reduces to [5] 

_dls sing d6 
M|rs ->°° " r| sin A^dAf 1 UJ 

Fig. [5] shows the magnification as a function of the the azimuthal shift [5] . 
All images have been put in sequence by adding (m — l)7r to the abscissa of 
images of order m. We see that the magnification diverges at all multiples 
of 7r. This divergence comes from the tangential magnification (|15p . which 
means that images of a finite size-source become large arcs and then merge 
into rings when the source is in front or behind the black hole. The primary 
image (0 < A<fi < it) has otherwise magnification very close to 1, whereas all 
other images have an exponentially decreasing magnification. The secondary 
image interpolates between the two behaviours. We conclude that the lines 
defined by A<p = and A<p = it and extending from the horizon to infinity 
(2M < r < oo) represent degenerate (with zero section) caustic tubes for the 
Schwarzschild metric. 



7 Here r)ij = diag(l, r 2 , r 2 sin 2 ■&) is the Euclidean metric in polar coordinates as the 
source is in the asymptotic region. 
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The parity of the images is determined by the factor sinZ\0 in the tan- 
gential magnification. All direct images have positive parity and all indirect 
images have negative parity. 

We conclude by saying that possible effects beyond the geometrical optics 
scheme for photons travelling close to the photon sphere are still awaiting a full 
investigation. The existence of nonlinear dispersion relations for the photons 
propagating close to the photon sphere could have important consequences 
not yet considered [16]. 



2.4 Analytical approximations 



The elliptic integrals of the exact treatment of light bending are not partic- 
ularly illuminating and prevent an analytical inversion of the lens equation 
(|12[) . For this reason, several approximation schemes have been developed in 
different limits. With their analytical formulae they shed new light on the plots 
of the previous subsection. Furthermore, taking advantage of the perturbative 
schemes in different limits, it is possible to compare different black hole met- 
rics coming from alternative gravity theories and use gravitational lensing to 
discriminate among different scenarios. 



2.4-1 Weak deflection limit 



The WDL amounts to setting w > M in the deflection integrals ([5]) or (fpQ)) 
and (fTl"]) . which also implies r m ,rs 3> M. We are thus in the lower-right 
corner of Fig. []] or the upper side of Fig. [3^. When A<f> is significantly less 
than 7r, the gravitational lensing effect on the primary image is negligible 
while the secondary image is not in the WDL regime. Therefore, the WDL 
approximation is only useful in the case A<p ~ 7r, where it can be nicely 
employed to describe the primary and secondary images. 

In the WDL it is possible to study a generic static spherically symmetric 
metric by using the PPN formalism. In isotropic coordinates, the metric can 
be parameterized as 



ds 2 = A'(r')dt 2 - B'(r') 



s 2 dn 2 



,M 



M 2 

AV) = 1-2^ + 2/3'^. 

r r' 

,M 3 „,M 2 



(21) 
(22) 
(23) 



In the Schwarzschild limit, the PPN coefficients become a' = f5' = 7' = 8' = 1. 

The integrals in the azimuthal shift (1101) can be easily expanded in powers 
of M. In particular, the deflection angle © has been calculated to the second 
order in M by Epstein and Shapiro already in 1980 [T7] (see also [Hj). Defining 

the Einstein angle 

0e = y/4Mr s /{r dos), (24) 
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a convenient expansion parameter is [19] 

e = 9 E d os /(4r s ). (25) 

Without gravitational lensing, the observer would see the source at angle /? 
from the lens. Simple geometry relates this angle to the azimuthal shift A<fr as 
dos sin (3 = rs sin Acj). Setting y = P/Oe, it turns out that to lowest order 9 ~ 
9e and we can expand the position of the images as 9 = 9 E (9 +e9i+e 2 92 + ■ ■ .). 

Expanding the exact lens equation (TT2l) , we can collect all terms according 
to their order in e. Finally, we can calculate the position of the images and 
their magnification /i = (/Xo + £A*i + • • •) order by order. We just quote here the 
final result^! found by Ebina et al. [3D] in the Schwarzchild case, by Sereno [21] 
for Reissner-Nordstrom and by Keeton and Petters for a generic PPN metric 

m 

9 = (y + vV + Ai) A e 1 = A 2 /{A l + iel) (26) 
fi = 166*/(169*-Al), fi 1 = -l69 3 Q A 2 /(A 1 + 49 2 ) 3 (27) 
with the coefficients A\ and A 2 fixed by the PPN metric as 

A x =2(c/+7'), A 2 = 2a' 2 - p + 2a V + 3<J'/4. (28) 

In principle, if the source is close enough to the lens so as to raise e signifi- 
cantly, by studying gravitational lensing we can gain access to the first order 
correction in the WDL expansion and thus infer the values of PPN parame- 
ters. Surprisingly, it turns out that all observables that are independent of the 
unknown source position are unchanged at first order in e [20] for all PPN met- 
rics agreeing with General Relativity at zero order (thus having a' + 7' = 2) 
[19] . This fact was prefigured by Damour and Esposito-Farese [35] using an 
independent field theory approach to gravity theories in which they showed 
that 2PN parameters do not enter the equation of motion of light. 

Finally, we note that the exact lens equation (TT21 is not really necessary to 
calculate the first order correction, as one would reach the same result using 
the classical small angles lens equation with a second order expansion of the 
deflection angle $Z§. The second order corrections in e, however, can only be 
calculated by the exact lens equation or by the Ohanian lens equation (TT^|) 23 , 
as other lens equations would deviate from the correct results. The Ohanian 
lens equation, instead, would fail at the fourth order only 13 . 

2.4-2 Strong Deflection Limit 

The expansion of the elliptic integral in Darwin's formula at its divergence 
at r = 3M leads to the SDL formula ©, which shows that the divergence at 
u = 3vo.M is logarithmic [TllTi rHll2~i] . This is not a peculiarity of Schwarzschild 
metric, as even in Reissner-Nordstrom metric a similar divergence can be found 



In these formulae y is positive for the direct image and negative for the indirect image. 
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[25] . Actually, it is possible to prove that the logarithmic divergence in the 
deflection angle is a universal feature of all static spherically symmetric metrics 
endowed with a photon sphere [26JH1- 

In order to calculate the SDL in a given metric, it is first necessary to 
calculate the radius of the photon sphere. This is obtained by imposing that 
the denominator of Eq. ^ has a double root at r — r m . This is equivalent to 
the equation 

C'(f)/C(f)=A'(f)/A(f), (29) 

which yields the radius of the photon sphere f [3] . The minimum impact pa- 
rameter u is then given by Eq. (|3J with r m = f. 

We can then expand the integrals in the azimuthal shift (fT0|) in terms of 
e = (u/u — 1) to get the SDL lens equation^] 

A4> + 2?itt = -a log[e/(l - r/r s )} + b, (30) 

where the coefficients a and b depend on the metric considered^ and on the 
source radial coordinate rg. Indeed, the SDL coefficients may vary in different 
alternative theories of gravity, representing a sort of identity card of the black 
hole metric and a possible key for discriminating between theories predicting 
deviations from General Relativity in strong fields 29 , 30 . In the Schwarzschild 
case 



1, b = -2 log 



(3 + Vtj (3 + ^3 + 18M/r s ) /(36V6) . (31) 



Corrections to Eq. Q30p at higher orders in e have been worked out only in the 
Schwarzschild case for r$ — > 00 [31J . 

This lens equation can be used to describe all images with n > 1, generated 
by photons looping at least once around the black hole (cfr. Fig. [3^). It becomes 
inaccurate only for the secondary image. It may even be used for sources inside 
the photon sphere (Fig. [3Jd) . The position of the images is found by inverting 
Eq. ((3"Tj) and remembering that 8 = u/ro in the limit ro S> M. We have [28] 

1 + (1 - f/rs)^-^- 2 "^/ 5 ! . (32) 

The indirect images are found by replacing A<p by 27r — Acj> as usual. This 
formula indicates that the higher order images tend to the shadow border 9 
with an exponential law (which justifies our choice of the logarithmic scale in 
Fig-©. The images form outside the border for r$ > f and inside it otherwise. 



9 It is easy to prove that all static spherically symmetric asymptotically flat metrics with 
an event horizon have a photon sphere. 

10 The SDL was developed in Ref. 1201 for standard lensing only. The extension to the retro- 
lensing configuration was done in Ref. 1271 . Here we present the most general treatment of 
Ref. [28] . 

11 They can be calculated analytically or numerically using the procedure described in Ref. 
26 for sources at infinity and in Ref. 28 for any source positions. 
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For sources very far from the black hole, the magnification is 5,26J 

M „ = ( £os\ 2 (S-^-a^/a (33) 

\ rs J asmA<p 

which nicely reproduces the behavior of Fig. [5] and also proves that higher 
order images are exponentially fainter and fainter. 

It is interesting to note that the time delay between two consecutive higher 
order images is approximately AT = 2it9ro, which could be used in principle 
to infer the mass and distance of the black hole [32] . 

Another interesting application of the SDL scheme is in the study of a static 
spherically symmetric black hole embedded in an external gravitational field. 
In the case of the standard WDL, Chang and Refsdal [33] showed that the zero 
order Einstein ring becomes elliptical and the corresponding degenerate caustic 
tube at A(p = tt acquires a diamond-shaped cross-section with angular size 
proportional to the external shear 7 and 9e- In the SDL, a correct description 
of the tidal forces 34 leads to the conclusion that also higher order Einstein 
rings gets slightly distorted and their corresponding caustic tubes are diamond- 
shaped with size proportional to 7$ [55] , 

2.4-3 Non-perturbative methods 

While the WDL describes the two main images in the best alignment regime 
between source and lens (A<f) ~ tt) and the SDL describes all higher order 
images (n > 0) , no perturbative approximations are possible for the secondary 
image when the source is not aligned behind the black hole (n = and < 
A<f> < tt). Indeed Figs. [3J and [5] show that the secondary image interpolates 
between the WDL and the SDL regimes when Acf> runs from tt to 0, going 

from e ~ e E to e ~ e. 

A nice strategy to obtain simple analytic formulae has been proposed by 
Amore and collaborators [36]. First, the integral ([2]) is written in the form 

l 

a = 2 f - dz - n (34) 

J y/V(l)-V(z) 

with z = r m /r and 

V(z) = z 2 C(r)/B(r) - C 2 (r)A(r m )/[B(r)A(r)C(r m )} + A(r m )/C(r m ) (35) 

being a sort of potential in the variable z. 

Then, the potential V(z) is replaced by a simpler form Vq and the integral 
becomes 

l 

r dz i _ ^ 

J y/V (l)-V Q (z) y/T+JAjz) 

w-vrcr^- 1 - (37) 

Vo(lJ - V {z) 



Interestingly, the lensing problem is reformulated as an oscillator problem. 
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Finally, the integral is solved as a series in powers of S and calculated in S = 1 
(where Eq. (j3"6")l coincides with (jM)) ). The convergence of the series depends 
on the chosen form of the auxiliary potential Vq- If Vo depends on some shape 
parameter A, this can be fixed by requiring the minimum sensitivity of the 
deflection angle (as in a variational method). 

As a practical example, taking Vq — \z 2 in the Schwarzschild metric, 
Amore and Arceo [36 obtained 

A/> = tt ( 1 = - 1 ) , (38) 

\y/l-8M/(nr m ) J' 

which works nicely in the intermediate regime for the zero order indirect im- 
age. More sophisticated potentials may lead to even better approximations 
englobing the main features of the WDL and SDL [33] at the price of dealing 
with more involved formulae. 



3 Rotating black holes 

Most astrophysical systems are known to have non-negligible angular momen- 
tum. Therefore it is natural to expect that physical black holes should be prop- 
erly described by the Kerr metric. The study of null geodesies in Kerr geometry 
has been initiated by Carter in 1968 [37], who separated the Hamilton- Jacobi 
equation and wrote the geodesies equations in terms of first integrals of the 
motion 

where 

= Q + a 2 cos 2 d — J 2 cot 2 0; A = r 2 - 2Mr + a 2 (41) 
R = r 4 + (a 2 — J 2 — Q)r 2 + 2M(Q + ( J - a) 2 )r - a 2 Q, (42) 

J and Q are constants of motion labelling the geodesic under examinational, a 
is the angular momentum of the black hole divided by its mass and expressed 
in natural units. The Kerr metric describes a rotating black hole if a < M with 
horizon at rji = M + \J M 2 — a 2 , otherwise it represents a naked singularity. 
Both f and $ may change sign during the motion. Therefore, every integral 
should be calculated from the emission point to the first inversion point (if 
any), then between each pair of consecutive inversion points, up to the final 
observation point. In the following, we will use the notation = cos F^l. 

13 J is the projection of the angular momentum of the photon on the rotation axis, Q is a 
combination of the square of the other components of the angular momentum of the photon 
with the black hole angular momentum. 

14 Not to be confused with the magnification of the previous section. 
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The main difference with respect to spherically symmetric black holes is 
that the motion of the photon is no longer confined on a single plane, as the 
angular momentum of the black hole induces an orbital precession around its 
rotation axis. As a consequence, the polar angle $ cannot be eliminated by 
a simple choice of coordinates and the problem is fully three-dimensional. In 
principle, all integrals can be written in terms of elliptic integrals (see e.g. 
[35]). This is not particularly illuminating for the analytic point of view, but 
may be useful for implementing numerical codes for Kerr-geodesics tracing [391 
I51?ll4"l"ll4"2] . In particular, the analytic deflection angle for photons confined in 
the equatorial plane takes a reasonably short expression |43j . 

A thorough presentation of the classification and the properties of Kerr 
geodesies is in Ref. [24 . In order to investigate the rotation of the polarization 
plane of electromagnetic radiation (not discussed in this text), it is useful to 
introduce a tetrad formalism |44) . As for spherically symmetric case, several 
studies of Kerr lensing in the WDL [17U55U5SU57] and in the SDL [2511551 
|4"5] have been carried out. These studies allow to get interesting analytical 
approximations in some limits that we will sometimes recall later. 

In the following subsections, we will give a sketch of the main characteristics 
of lensing by rotating black holes. We will first describe the shadow of the black 
hole, then we will discuss the caustics and finally the formation of the images. 



3.1 Shadow of the Kerr black hole 

Null geodesies are fully identified by the constants of motion J and Q. For a 
static observer in the asymptotic flat space (ro 3> M), the position in the sky 
in which a photon is detected is completely determined by the values of these 
constants. In fact, defining fip = cos$o, we have 



where Q\ is the angular distance from the black hole center projected orthogo- 
nally to the spin direction and 8 2 is the same angular distance projected along 
the projection on the spin on the sky. The problem of determining the position 
of the images of a given source can thus be reformulated in terms of J and Q. 

Similarly to the spherically symmetric case, the integrals in the radial coor- 
dinate diverge logarithmically if the argument of the square root has a double 
root. This happens when R and dR/dr simultaneously vanish at the same 
point f. The two equations R = and R' = can be solved in terms of the 
two constants of motion J and Q as functions of f, so as to find those geodesies 




(43) 



62 = ±To 



(44) 
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Fig. 6 The shadow border as seen by an equatorial observer for different values of the black 
hole spin. 6\ and 02 are the angular coordinates in the observer's sky in units of 2M/ro- 
The black hole is in the origin and the spin axis points in the positive 02 direction. 



which asymptotically approach an orbit at a fixed radius r. Explicitly, we have 

_ f_Hr-3M)+a^M + r) 

f 3 \l6M 3 a 2 - f(f ~ 3M) 2 ] 
Q & = L 16MMf-M)* (46) 

which describe a locus in the (J, Q) space parameterized by f . Only positive 
values of Q describe geodesies approaching infinity (see Eq. (04])). This con- 
strains the possible values of f into a finite interval [f + ,f_] defined by the 
conditional Q > 0. The two extreme values f± lead to a vanishing Q, which 
corresponds to a photon moving on an equatorial orbit. The lower value f+ 
corresponds to a photon moving with positive J, i.e. rotating in the same 
sense with the black hole (prograde photon) . The higher value r_ gives a neg- 
ative J, thus describing a photon counter-rotating with respect to the black 
hole (retrograde photon). In-between, we can have orbits with any possible 
inclinations. 

From the point of view of the observer at infinity, thanks to Eqs. (|43[) 
and (|44p . all geodesies asymptotically starting in an orbit at fixed radius f, 
described by Eqs. (|4"5|) and P5|) . also describe a locus in the sky coordinates 
(Ox, 82). Notice that if the observer is not on the equator, equatorial orbits 
will never reach it. Indeed, the range of values of f that describe orbits ending 
at the observer is defined by the condition that the argument of the square 
root in Eq. (|44[) is positive. The new range [f<,r>] is obviously smaller than 
[r+, f_] and coincides with it if the observer is on the equator. 



The third degree equation Q = has always two roots outside the horizon. 
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Analogously to spherically symmetric black holes, the locus (8~i(f), 62(f)) 
is called the shadow border of the Kerr black hole. In Fig. [5] we show the 
shadow border for three values of the spin parameter a. Indeed, when a — > 0, 
the shadow becomes a circle with radius 3*/3M/ro as for the Schwarzschild 
black hole. When we turn the spin on, the shadow shifts to the right (the spin 
direction points to the top in Fig. [5]) and becomes more compressed. With our 
coordinate system, photons detected on the left side of the black hole have 
rotated in the same sense of the black hole whereas photons detected on the 
right side have counter-rotated. As mentioned before, prograde photons are 
allowed to get closer to the black hole (up to r+) and therefore reach the 
observer with a small angular separation from the black hole. In the extremal 
limit a — > M, the shadow border on the left (prograde) side is flattened while 
f+ equals the horizon radius rjj. Beyond this limit, the horizon disappears 
and the Kerr metric describes a naked singularity whose appearance is no 
longer determined by the last circular orbits of the photons but depends on 
the physical interpretation of the singularity itself 50 . 

The idea of measuring the spin and other characteristics of astrophysical 
black holes from the shape of the shadow border dates back to Bardeen in 
1972 51j and has been re-proposed several times [24,50,52 . However, the real 
appearance of an astrophysical black hole is largely dictated by the properties 
of the accretion flow, as we shall discuss in the next section. 

3.2 Caustics 

As customary in gravitational lensing, the general functioning of a lens model 
is determined by its caustic structure. In fact, by studying caustics, we can find 
the points in which a source is maximally amplified and the regions in which 
new pairs of images are created. As showed in the previous section, spherically 
symmetric black holes have degenerate caustics in the form of zero-size tubes 
lying on the observer-lens axis behind and in front of the black hole. It can be 
easily imagined that the loss of spherical symmetry makes these tubes shift 
and acquire a finite extension. The existence of finite-size caustics in the Kerr 
metric was pointed in 1972 by Cunningham and Bardeen |53j . who studied 
gravitational lensing of a star orbiting a Kerr black hole. They found that 
additional pairs of images were created when the source was in some particular 
regions. The complete shape of the primary caustic was instead illustrated for 
the first time by Rauch and Blandford [38] who showed that it becomes a 
finite-size tube with astroid cross-section, winding around the horizon. Lately, 
we discussed the complete caustic structure of the Kerr metric, illustrating 
higher order caustics and their dependence on spin and the observer position 

In order to calculate caustics, it is necessary to put Eqs. (j39|) and (|40|) in 
the form of a lens mapping relating the source coordinated (fisi&s) to the 



16 The source radial coordinate rg is considered as a parameter. For each value of rg we 
then obtain a section of the caustic tube at fixed distance from the black hole. 
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image position (61, 62)- After that, we have to calculate the Jacobian determi- 
nant and find the locus defined by its vanishing (critical curve). Finally, the 
critical curve is mapped to the source coordinates through the lens mapping 
to obtain the caustic. Everything can be done numerically, but the details are 
quite cumbersome and are summarized very shortly in this treatment. On the 
other hand, some analytic formulae can be obtained in the WDL and SDL 
analogously to what happens in spherically symmetric black holes. 

In brief, it is convenient to replace the coordinates in the observer's sky 
(01,02) by pseudo-polar coordinates centered on the black hole position. Fix- 
ing the angular coordinate 77, the angular distance e from the shadow, defined 
in analogy to what is done in the SDL, determines how much the photon is 
deflected from the black hole. In fact, a photon detected far from the black hole 
(e ^> 1) is weakly deflected, whereas the deflection diverges when the photon 
is detected by the observer infinitely close to the shadow (e — > 0). The number 
of oscillations in the polar angle i? depends on the radial coordinate e. It turns 
out that in the range [e m ,e m +i] in which the photon performs m inversions 
in the polar motion, there is only one point in which the Jacobian determi- 
nant vanishes. Therefore, we can identify the caustic order with the number 
of inversions in the polar motion of photons that form the corresponding crit- 
ical images. The primary caustic is formed by photons experiencing only one 
inversion in the polar motion. The secondary caustic is formed by photons 
performing two inversions in the polar motion and so on. 

Fig. [7^ is a 3-dimensional representation of the primary caustic |38] for a 
quasi-extremal Kerr black hole and Fig. [7j3 shows its cross-section at various 
distances. We can see that the zero-size tube of the Schwarzschild black hole 
here becomes a tube with astroid cross-section winding clockwise around the 
horizon an infinite number of times [38(154] ■ At large distances (rs 3> M), the 
caustic can be described analytically in the WDL, since the photons of the 
critical images are weakly deflected. As shown by Sereno and De Luca [47] . 
the primary caustic can be written in the coordinates (/zg = cos?9s, (f>s) as 



with a m = a/M, rj ranging from — tt to tt and e being the WDL expansion 
parameter already introduced in spherically symmetric black holes in Eq. (|25[) . 
It is interesting to note that the caustic remains infinitely thin at second 
and third order in e. We need to push the expansion to the fourth order in 
order to find a finite extension with the typical astroid cross-section. This fact 
signals that to lowest order the Kerr black hole can be perfectly mimicked by 
a Schwarzschild black hole with a finite shift, as previously noted by Asada 
and Kasai ,46]. The angular size of the primary caustic varies with distance 
as e 4 ~ r^ 2 . Therefore its physical extension shrinks as r^ 1 far from the black 




(t>s = -7T - 4a m e 2 - - 





(47) 



(48) 
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Fig. 7 (a) The primary caustic tube for a = 0.9998M and jio = in a 3-dimcnsional repre- 
sentation, with pseudo-euclidean coordinates x = r sin cos <f>, y = r sin $ sin <j>, z = r cos i?. 
The spin axis is directed toward the top. The straight solid line indicates the direction 
towards the observer, whereas the dashed line points in the opposite direction. The pri- 
mary caustic for a Schwarzschild black hole (a = 0) coincides with this dashed line, (b) 
Cross section of the same primary caustic at various distances. From left to right, the 
radial coordinate is r = 5M, 4.5M, 4M, 3M, 2.3M, 1.82M, 1.42M, 1.22M. (c) 3-dimensional 
picture of the second order caustic surface of an extremal Kerr black hole for an observer 
on the equatorial plane in the direction indicated by the solid line. The surface has been 
plotted for radial distances in the range [2.2M, 20M]. (d) Second order asymptotic caus- 
tic for an equatorial observer /i = and different values of the spin. From left to right, 
a = 0.02M, 0.2M, 0AM, 0.6M, 0.8M, 0.9998M. 

hole, severely limiting chances of observations of spin effects in the WDL. 
The effects of the spin becomes sizeable at stronger deflections, which can be 
obtained either with sources at distances from the black hole comparable with 
the Schwarzschild radius or considering higher order images. 

The secondary caustic involves photons with two inversions in the polar 
motion. In the Schwarzschild limit it reduces to an infinitely thin tube on the 
same side of the observer. Its corresponding critical curve is the first retro- 
lensing Einstein ring. In Fig. [7J1 we can see what happens when we turn the 
spin on |54j . The secondary caustic becomes very large and for quasi-extremal 
spins has an angular extension larger than tt in the azimuthal direction. Fur- 
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thermore, the angular size does not decrease with the distance rg, and the 
whole cross-section tends to an asymptotic astroid shape that is what we 
show in [Jji. An impressive 3-dimensional view of the secondary caustic in the 
quasi-extremal case is shown in Fig. [7b. The caustic surface includes most of 
the equatorial plane assuming a discoidal shape. Creation and destruction of 
images is thus much more common in the second Einstein ring, than in the 
primary one. The spin effects, which are crucial to determine the size and the 
position of the secondary caustic, are then much more evident on higher order 
images, but we must be aware that these are much fainter than the primary 
images^]. Also for higher order caustics it is possible to derive analytical ap- 
proximations in a specific limit (SDL and small values of the spin, with rs free 
to vary from r# to infinity) [2"8"ll4"9"] . In this case for a caustic of order m > 2 
we have 



MS = (-I)™ [/iO + Rrn(l - Mo 2 ) 3/2 Sin 3 T, 
4> S = -miT - A(f> m + R m \/l - Ho 2 COS 3 7], 

where the azimuthal shift is 



A(/) m = -a ■ 



2miT 
3\/3 



2 log (2V3-3) +log 



(2y^+ V6Af + r s ) 
3y/(r s - 2M) 



(49) 
(50) 

, (51) 
• (52) 



As for the primary caustic, the shift is linear and the size is quadratic in the 
spin. It should be noted that for small values of the spin there is a degeneracjfjf] 
between the spin a and the latitude of the observer on the equatorial plane. 
The caustics become smaller and smaller as the observer is moved from the 
equatorial plane to the polar axis, finally returning to zero shift and size when 
the observer is perfectly aligned with the rotation axis restoring the axial 
symmetry. The degeneracy between spin and observer latitude is only broken 
at high values of the spin [54] . 



and the semi-amplitude of the caustic is 

Rm = i 



1 (5™ + 8V3 - 36) + (^ + 2r s -2^V6MT^) 
18 3\Z*yrsV6M + r s 



3.3 Formation of the images 

The discussion of the images in Kerr black hole lensing may start from the 
Schwarzschild limit. In this limit, we know that we have an infinite sequence of 
direct images and an infinite sequence of indirect images. We have noticed that 

17 Fig. [5] in the Schwarzschild case may be considered as a good reference even in the Kerr 
case for a rough estimate of the luminosity of images of various orders. 

18 The degeneracy is evident when switching to coordinates centered on the zero-order 
caustic position fi = (—l) m fio- Then all quantities depend on the combination au 1 — 
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the image of order to has an azimuthal shift A<p in the range [(m — l)7r, mir]. 
In the Kerr case, it is possible to maintain this identification using the number 
of oscillations in the polar angle, instead [53] . 
Solving Eq. (f3T))) explicitly in /is, we have 



where k and tp are functions of J and Q, -K"(fc) = F(n/2, k) is the complete 
elliptic integral of the first kind and sn is the Jacobi elliptic function. As sn 
has period 4K(k), ip counts the number of half-oscillations on the equatorial 
plane. It turns out that there is always at least one image in each interval 
m — 1 < tp < m. More specifically, when the source is inside the caustic of 
order to, there are three images, otherwise only one. 

Comparing with the Schwarzschild black hole, in that case we had an image 
in each interval (m — l)ir < A<f> < rmr. Direct images were those with odd 
to and indirect images were those with even to. In the Kerr black hole the 
situation is similar, but the azimuthal shift is no longer suitable to label them, 
because it evolves differently for photons with different inclinations on the 
equatorial plane. The argument tp of the sn function in Eq. (1531) provides the 
correct variable. Its explicit expression is reported in Ref. |54j . 

In Fig. [8] we show the position of the images of order 3 and 4 for a source 
very close to the caustic of order 3, which involves photons performing a bit 
more than one loop and a half (\<fi\ > 3ir). When the source is outside the 
caustic (Fig [5^) , there is one third order image (on the left) and one fourth 
order image (on the bottom right). The images appear as very stretched thin 
arcs. When the source is inside the caustic, two more third order images appear 
on the top right, while the fourth order image remains lonely. 

The Kerr lens equation for sources close to caustics assumes a very simple 
form for small values of the spin. In this limit, we have [3Ulf4^1l4T] 



where (W SI S(f> s ) is the source position relative to the center of the caustic, 
R m is the caustic radiuJ^l. and rj is the angular coordinate around the shadow 
border in the observer's sky as defined in Section 13.21 The distance from the 
shadow e remains the same as in Schwarzschild to lowest order in a. 

For near-extremal black holes (a — > M) , the azimuthal size of higher order 
caustics exceeds 2ir [53]. As <p is a periodic coordinate, this means that we 
should consider the possibility that a source crosses the same caustic surface at 
different number of loops. When this happens, new pairs of additional images 
appear at the same order to but with A(j) differing by multiples of 27r. From 
the physical point of view, this means that these pairs of images are formed by 
photons performing the same number of oscillations m in the polar angle but 
a different number of loops in azimuth 4>. These new pairs of images appear 
on the flattened side on the left of the shadow and are formed by prograde 
photons moving very close to the equatorial plane. As the extension of higher 



/is = ±/i + sn (2K{k)ip, k) , 



(53) 



Sd s cosr] + (— I) m (<5</> s — R m cos ij) sin?? = 0, 



(54) 



It can be obtained from Eqs. (1476 and 1148 II for m = 1 and from Eq. H52II for m > 1. 
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Fig. 8 Images of order 3 and 4 of a source of radius 2M at rg = 200M seen by an observer 
on the equatorial plane of a black hole with spin a = 0.6M. In panel (c) we show two 
positions of the source with respect to the third order caustic. In panel (a) we show the 
images corresponding to the empty circle along with the shadow border and in panel (b) we 
show the images corresponding to the filled circle. 



order caustics increases exponentially with the order m, the number of images 
also grows exponentially with m in extremal Kerr black holes. However, the 
magnification drops even faster, so that it will be extremely difficult to observe 
this effect. 



4 Observational perspectives 

Black holes offer a possibility (unique in the astrophysical panorama) of study- 
ing light bending beyond the first order WDL approximation, which governs 
all known cases of gravitational lensing. However, unveiling this new phe- 
nomenology requires very powerful technical capability. Let us start with a 
general discussion. 
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The first corrections to the primary and secondary images, which are well- 
known in all standard gravitational lensing situations, scale with the combi- 
nation 9e£, which is proportional to M/ro- Interestingly, also the angular size 
of the shadow of the black hole and all observables related to higher order im- 
ages are proportional to the ratio M/ro- We shall therefore take the shadow 
angular size 9 = 3y/3~ M/ro as the guiding quantity to evaluate astrophysical 
black holes as candidates for gravitational lensing beyond the WDL standard 
approximation. 

Typical stellar black holes in the Galaxy have 

Osteiiar = 5.1 x 10~ n arcsec \ (i~£p~c") ' *' 55 ' ) 

which is definitely a too small angle to be observed. Much promising candidates 
are massive black holes at the center of galaxies. In particular, the black hole 
in the center of the Milky Way, identified with the radio source Sgr A* [35] , has 
mass 4.31 xl0 6 M Q [SB] and lies at about 8 kpc from the sun. The angular size 
of its shadow is thus 0s gr A* = 28yuas. It is worth noting that some supermassive 
black holes in other galaxies reach similar orders of magnitude. For example, 
M87 has a central black hole of 3.2 x 1O 9 M at a distance of 15 Mpc from us 
[57] . resulting in a shadow radius of ll^tas. However, we must be aware that 
the luminosity of the sources falls down as d,Q 2 s , limiting possible gravitational 
lensing observations. The black hole in the center of our Galaxy was proposed 
as the best candidate for black hole lensing by Cunningham and Bardeen [53] . 
Most of the works in this field have thereafter focussed on this nearby massive 
black hole. 

Sgr A* has been intensively studied in the radio and X bands, showing the 
existence of daily flares coming from the inner accretion regions and probably 
generated by clumps of matter falling down the horizon. As the Galactic center 
region is optically thick, the near infrared band offers the best opportunity to 
study stellar sources around Sgr A*. The latter has only been detected in its 
flare state in the K-band at 15 magnitudes [58]. We refer the interested reader 
to the textbook by Melia [59] for an extensive review summarizing observations 
and models of Sgr A* and its environment. 

Weak gravitational lensing by Sgr A* on background stars should be quite 
common [50], though it has never been identified because of the crowded en- 
vironment. Unfortunately, spin effects are negligible on the primary and sec- 
ondary images in the WDL regime. This can be deduced by the size of the 
primary caustic, which shrinks to zero at large distances. In the case of Sgr A*, 
the size of the primary caustic is only 1170 km at r$ = lOORsch- It is therefore 
necessary to consider either sources close to the black hole or higher order im- 
ages, which are however very faint and appear very close to the shadow of the 
black hole within its accretion disk [51161). In the following subsections, we will 
discuss several observational proposals which aim at detecting effects beyond 
the basic WDL in gravitational lensing by Sgr A* or other black holes. 
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Fig. 9 Apparent positions of the primary (solid line) and secondary (dashed line) images 
of a source at r$ = 20Af from an extremal Kerr black hole. The observer is at a polar angle 
•So = 84.24°. Ticks mark the positions of the images at 10 equally spaced times. As the 
source crosses the secondary caustic, two new second order images are formed at the point 
© and disappear at (gi when the source exits from the caustic. The inner dashed circle has 
angular radius M/ro- Picture taken from Cunningham and Bardeen (1973) |53l . 

4.1 Sources orbiting the black hole 

Stars orbiting the black hole provide natural sources for gravitational lensing. 
This possibility was first considered by Cunningham and Bardeen |53j . who 
traced the apparent trajectories in the sky of the primary and secondary images 
for a star orbiting a black hole in a circular orbit. In Fig. [9] we see that the 
primary image is practically unaffected when the source passes in front of the 
black hole. When the source transits behind the black hole, the primary image 
moves around the primary Einstein ring, designing a hat figure. Meanwhile, the 
secondary image is opposite to the primary image slightly inside the primary 
Einstein ring. However, when the source moves in front of the black hole, the 
secondary image is closer to the black hole, approaching the secondary critical 
curve. Very interestingly, this picture shows the formation of two additional 
second order images and the annihilation of one of these with the original 
secondary image. These two events occur when the source enters and exits the 
giant secondary caustic showed in Fig. [7b. 

We do not know of any stars on a circular equatorial orbit so close to Sgr 
A*. Yet, the picture of Fig. |H] may be likely applied to a bright spot of hot 
material lying on the accretion disk. The existence of such features has been 
proposed by Abramovicz et al. in order to model the rapid X-ray variability 
of Active Galactic Nuclei (AGN) [62]. Substantial amount of work on this 
possibility has been developed thereafter 38 , 39 , 63 . Assuming that we do not 
have enough angular resolution to see the distinct images, the flux of the hot 
spot would be periodically modulated through gravitational lensing by the 
black hole. As the spot would rotate very quickly in a strong gravitational 
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Fig. 10 Light curve of a bright spot orbiting a quasi-extremal Kerr black hole at different 
radii as indicated in the curve labels. Time is normalized to the period. Picture taken from 
Karas, Vokrouhlicky and Polnarev (1992) I39| . 



field environment, Doppler and gravitational red-shift cannot be neglected in 
a correct treatment. Note that if the disk is opaque to the electromagnetic 
radiation considered, only the primary image contributes to the observed flux, 
whereas the second order and all higher order images are stopped by the disk. 
In Fig. I10[ we see that the light curve is dominated by a peak occurring when 
the spot is close to the primary caustic behind the black hole. This peak may 
be followed by a secondary peak, depending on the spot distance, caused by 
Doppler enhancement. 

Gravitational lensing of hot spots is very involved because of the large num- 
ber of parameters to be guessed, such as the size, the shape, the position and 
the internal velocities of the spot, besides the disk influence on the electromag- 
netic radiation emitted by the spot. On the other hand, we know of many stars 
orbiting Sgr A* at larger distances on eccentric orbits [56,64 . Spectra of these 
stars show that they are early-type unexpectedly young stars; thus, their ori- 
gin is still under debate. If lensed, these stars would provide much cleaner data 
|65j . The best candidate is S6, whose secondary image reaches K — 20.8, with 
an angular separation of 0.3 mas from the central black hole [66]. This is just 
beyond the limits of the forthcoming interferometric instrument GRAVITY at 
VLT [57]. 

Another interesting possibility comes from Low Mass X-Ray Binaries (LMXB) 
orbiting Sgr A* [55]. In these objects, the X-ray radiation is emitted from a 
very small region, providing good candidates for probing the caustic structure 
of the black hole in the case it has a non- negligible spin [49] , 



26 



4.2 Broad Iron lines in the X-ray spectrum 

Among the known manifestations of gravitational lensing by black holes, the 
broadening of the iron lines is the one that has attracted more interest, as it 
is accessible with present technology. This phenomenon was first discovered 
in the Seyfert 1 galaxy MCG-6-30-15 by Tanaka et al. [59], in which a bright 
emission line is present over the X continuum at 6.4 keV. This line has a 
characteristic doubly horned structure, with the blue horn higher than the red 
horn. 

This shape has been interpreted as follows: cold material in the accretion 
disk is ionized by radiation from hot material surrounding it. The following 
recombination is responsible for the emission line, whose broad structure is 
due to the strong rotation of the disc, which Doppler shifts the emission. The 
two horns come from material whose velocity is directed along (or closely to) 
the line of sight in both ways, with the blue horn appearing more intense 
thanks to Doppler boosting. This picture also applies to several lines at lower 
energies, which are generated far from the black hole. The specifity of the iron 
lines comes from the fact that they probe very high energies. In particular, 
the K a line at 6.4 keV comes from electrons falling directly to the K-shell. 
Such energies are only reached in regions very close to the black hole (few 
Schwarzschild radii). At such distances, general relativistic effects cannot be 
neglected in the description of the line profiles and the whole spectrum in 
general [70] , 

The iron K a line has been observed in several AGNs and the number of 
works attempting a numerical or analytic description of relativistic iron line 
profiles has flourished. The first treatments in the Schwarzschild metric have 
been done in Refs. [71]. Lines in the Kerr metric have been studied in Refs. |72) . 
Usually, only the primary image of the disk fluorescence is taken into account, 
because the accretion disk is supposed to be thick. However, the second order 
image is typically visible if the disc is seen nearly edge-on. The relevance of 
higher order images has been pointed out in Ref. [73] and later in Ref. [IT] . 
A non-exhaustive reference list for later studies is in Ref. [741 . Recently, also 
the L-line at 0.9 keV, coming from electrons absorbed to the L-shell, has been 
observed [75] , 

Fig. ITT1 shows some relativistic line profiles calculated by Fanton et al. [41] . 
The lines are broader for a disc seen edge-on, showing up to four subpeaks. 
The line returns single-peaked and is slightly redshifted for a face-on disc. The 
lines are broader in a Kerr black hole, but there is almost no difference if the 
inner radius of the emitting region is larger than 6A/. 

Models of relativistic line profiles depend on many parameters and as- 
sumptions on the underlying model, such as opacity, thickness, homogeneity, 
outer and inner radius of the emitting region, emissivity profile and more and 
more. It is therefore difficult to use them to extract the black hole parameters 
with confidence. Nevertheless, they witness the existence of accretion disks at 
few Schwarzschild radii of distance from the black hole and, at present, they 
represent the cleaner information we can get from these mysterious regions. 
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Fig. 11 Relativistic line profiles for emitting disks around a Schwarzschild (left panels) and 
a quasi-extremal Kerr (right panels) black hole. In the top panels, line profiles for various 
viewing angles <Dq are shown. In the bottom panels, the outer radius of the disk r ou t is 
varied. Picture taken from Fanton et al. (1997) |41| . 



4.3 Direct imaging 

Direct imaging of a black hole with a resolution of the order of the angular 
size of the Schwarzschild radius is one of the ancient dreams of astrophysicists. 
Up to now, we have always dealt with indirect information on the regions sur- 
rounding the horizon. For example, the measured spectra are the superposition 
of signals coming from different regions experiencing different physics, from the 
interstellar environment to the accretion flow up to possible jets. Direct imag- 
ing would disentangle these signals and allow a separate study of the different 
components, unveiling the secrets of the black hole environment. 

The first studies on the realistic appearance of a Schwarzschild black hole 
surrounded by a Keplerian disk date back to Luminet [14] . who showed how 
the disk on the prograde side of the black hole appears brighter thanks to the 
Doppler boost (which was also responsible of the asymmetry in the iron line 
profiles, as we have seen in the previous subsection). Fukue and Yokoyama 
drew pictures at different wavelengths |76) . Viergutz considered a vertically 
extended disk around a Kerr black hole [40], whereas the first pictures of 
the realistic appearance of a Kerr black hole were presented by Fanton et 
al. [H]. Afterwards, Falcke, Melia and Agol presented pictures of spherical 
inflow models taking into account the electron scattering in the mm band [77) . 
Polarimetric images were presented in Ref. [78] . A non-exhaustive list of later 
references is in Ref. |79j . 
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Fig. 12 Top panels: extremal Kerr black hole with radial inflow. Bottom panels: 
Schwarzschild black hole with flow on Keplerian shells. Left panels: unscattered plots. Mid- 
dle panels: the same plots including scattering by interstellar medium at A = 0.6 mm. Right 
panels: the same plots at A = 1.3 mm. Picture taken from Falcke, Melia and Agol (1999) 
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Iri Fig. [T^] we report the pictures by Falcke, Melia and Agol [77]. The dif- 
ference between Kerr and Schwarzschild is mostly in the shape of the shadow, 
which is flattened on the prograde side for Kerr. The prograde side of the 
flow is Doppler-boosted in both cases, giving poor indication on the spin of 
the black hole itself. We can also appreciate the blurring caused by scattering 
on electrons in the interstellar medium. Very Long Baseline Interferometry 
(VLBI) observations in the cm band show that the structure of Sgr A* is 
washed out by scattering, with its apparent size scaling with the wavelength 
as A 2 [SO]- Its intrinsic structure should be accessible at sub-mm wavelengths, 
on which Ref. [77] focussed. 

The best angular resolution so-far achieved by human instrumentations 
has been reached at A = 3 mm by VLBI on an intercontinental scale, with 18 
/xas [5_T]. Even if this is less than the apparent shadow of Sgr A*, interstellar 
scattering is still too high at this wavelength. Observations of M87 at this 
wavelength have also been carried out [55]. This has inspired the idea that 
the black hole in M87 could be an even more interesting target than Sgr 
A* [83], thanks to its jet structure. Meanwhile, projects for pushing the global 
network of VLBI to shorter wavelengths are under discussion [51] and represent 
the most promising frontier for direct imaging of black holes. Indeed, present 
observations at A = 1.3 mm seem to point that the intrinsic size of the emitting 
region of Sgr A* is smaller than the shadow size [53] . This is compatible with 
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a black hole surrounded by an accretion disk whose prograde side luminosity 
is enhanced by Doppler boost (see Fig. [T3]) 86 . 

Dynamical features such as hot spots rotating around the black hole may 
be directly imaged in the future [57]. These could also be detected at worst 
resolutions by following the motion of the centroid of Sgr A* [55] . 

In the near infrared bands, the best resolution is now obtained by the 
Keck interferometer and the VLTI. The latter will be upgraded with the new 
instrument GRAVITY [67] . which will take the resolution to 3 mas in the ri- 
band. This is far from the shadow scale of the black hole but is a great advance 
for the determination of stellar orbits, which put the tighter constraint to the 
mass of Sgr A*. It is also important for the study of flares in the infrared. 

Finally, the MAXIM project for developing X-ray interferometry on space- 
craft has been under study for several year^|. The incredible difficulties posed 
by this challenge will be rewarded by a black hole imager with resolution even 
better than 1 [i&s. The new frontiers opened by such a futuristic instrument 
can hardly be imagined. 
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